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Abstract 

A group symmetry analysis of the low lying levels of the spin- 1/2 kagome 
Heisenberg antiferromagnet is performed for small samples up to N = 27. 
This new approach allows to follow the effect of quantum fluctuations when 
the sample size increases. The results contradict the scenario of "order by 
disorder" which has been advanced on the basis of large S calculations. A 
large enough second neighbor ferromagnetic exchange coupling is needed to 
stabilize the \/3 x y/3 pattern: the finite size analysis indicates a quantum 
critical transition at a non zero coupling. 
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I. INTRODUCTION 



There are very few simple two dimensional magnets which fail to order at T = 0. There 
is now a large amount of evidence that the 5=1/2 nearest neighbor Heisenberg antifer- 
romagnet is ordered not only on the square lattice |l| but also on the triangular lattice 
(THA) The reduction of the order parameter by quantum fluctuations is about 40% 

for the square lattice and of the order of 50% for the triangular lattice: frustration indeed 
enhances the role of the quantum fluctuations but the relatively high coordination number 
plays against them. The kagome lattice which can be seen as a diluted triangular lattice (see 
Fig. p]) exhibits both frustration and low coordinance number and the S = 1/2 Heisenberg 
antiferromagnet on the kagome (KAH) is a good candidate for a disordered two dimensional 
quantum liquid. 

Exact diagonalizations on periodic samples have shown that the spin-spin correlations 
decrease indeed very rapidly with distance [|1C| -|T2]1 • Series expansion from the Ising limit and 



high temperature series point to the absence of magnetic order [|T^,|Tj]. Large-N approaches, 
using Sp(N) bosons fi~5| , find a disordered ground-state with unbroken symmetry for small 
enough spin while, using SU(N) fermions, a spin-Peierls phase or a chiral phase is suggested 
fT6|j . The best variational energies for the N=36 sample are built from the short-ranged 
dimerized basis 11171 . 



On the other hand, semiclassical approaches [|l8H 22"ll plead in favor of a magnetic order 
(the \/3xV3 state, see Fig. P induced by quantum fluctuations ("order by disorder" [f23|D . 
This kind of phenomenon has already been seen in much simpler situations where there 
is a classical degeneracy between two kinds of orders f24]j25l| . The Ji — J2 model on the 
triangular lattice for large enough J 2 (J2/J1 > 1/8) possesses classically degenerate ordered 
ground-states with respectively two and four sublattices | 2~6[| . Large spin calculations have 



predicted that quantum fluctuations should select the two sublattices order p7|j28| . The 
study of small samples spectra has confirmed this prediction and shown the mechanism of 
this selection 
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The kagome antiferromagnet is a more problematic issue because of the infinite num- 
ber of classically degenerate ground-states [[I(],|31[]. The linear spin wave approach for the 
kagome antiferromagnet does not lift the extensive degeneracy of the classical ground-state 
(at this order the spectrum of magnetic excitations possesses a whole branch of zero modes 
[P~5| , |5D|J5T| ] ) . One has to invoke non linear processes to stabilize the a/3 x a/3 ordered solution 
li,[20|,|2i. On the other hand, one should note that in this class of problem there exist some 



situations where the "order by disorder" phenomenon fails for systems with a ground-state 
manifold associated with an extensive entropy. An example is the quantum Heisenberg anti- 



ferromagnet on the triangular Husimi cactus |32| where the system prefers to remain a spin 
liquid rather than localizing in a particular ordered state that breaks the degeneracy of the 
classical ground state manifold. Moreover, the possibility of a quantum tunnelling between 
different classical states might prevent the system from localizing in a particular ordered 
solution [0,0 . 



This paper is a scrutiny of the conjecture of selection of order by quantum fluctuations 
based on the study of the low lying levels of exact spectra. The paper is organized as 
follows: A brief presentation of our method to investigate the finite size properties of an 
ordered antiferromagnet is given in the next section. In Sec. Ill, we apply this method to 
the search of the selection of the \/3 x a/3 state by quantum fluctuations. In Sec. IV, we 
introduce a second neighbor ferromagnetic exchange coupling (J 2 < 0) and show that for 
large enough | J 2 |, the system orders through the \/3 x a/3 state. The rigidity of the KHA 
and the hypothesis of an incommensurable planar order are investigated in Sec. V. Section 
VI, finally, summarizes our results and lists some open questions. 



II. REVIEW OF THE METHOD 

The presence of order in a quantum antiferromagnet is readily seen by examination of the 
low lying levels of its spectrum. The method used in this work has been described in details in 
^rjJ2^] , we give here a summary of its most important features. The first characteristic of an 
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ordered antiferromagnet (classical or quantum) is the existence of ferromagnetic sublattices: 
the total spin of each ferromagnetic sublattice is the collective variable relevant for the 
description of the low energy spectrum of the system. On a finite lattice with p sublattices 
of Q sites (the total number of sites being N = pQ), the p identical spins Q/2 couple to 
form the rotationally invariant states that are the low lying eigenstates of the hamiltonian. 
These p coupled spins form a "generalized top" which precesses freely: elementary mechanics 
(confirmed by the non- linear sigma model approach [p5l~|37D indicates that the leading term 
of this precession is: 

He, = (1) 
Z1 N 

In an ordered system, In, the inertia of the top is an extensive quantity, proportional to 
the perpendicular susceptiblity of the magnet P,|37|. In a disordered system with a gap 



this quantity is asymptotically constant, and at the quantum critical point between this 
two regimes In is expected to vary as iV -1 / 2 [S7|. The coupling of the pQ/2 spins gives 



Ns(N, S,p) states for each S value of the total spin of the sample and a total of (Q + l) p 
levels obeying the effective dynamics of Eq. ([!]). These levels form the quantum counterpart 
of the classical ground-state: they represent the tower of states indicated by Anderson in 
his 1952 famous paper In our spectral representation where the eigenstate energies are 
displayed as a function of S(S+1) in order to exhibit the free precession (Eq. ([!])), the tower 
is in fact a "Pisa Tower" with a slope decreasing with the sample size. In the thermodynamic 
limit, all these low lying levels (called QDJS in Ref. [|J for quasi-degenerate joint states) 
collapse to the absolute ground-state as iV -1 . An experimental proof of these levels has 
recently been reported in the analysis of macroscopic quantum coherence in antiferromagnets 
||39|| . It could be noticed that in the thermodynamic limit these levels give a ground-state 
entropy proportional to ln(iV p ). Above these first family of levels, the spectrum of an 
ordered antiferromagnet exhibits the magnon spectrum; the low lying levels of this part of 
the spectrum collapse more slowly than the QDJS to the ground-state with a scaling law 
jy-1/2 Three features of the "Pisa Tower" of QDJS are thus essential: 



i) the overall effective dynamics of a finite family of levels (of the order of N p ) and its finite 
size scaling leading to a clear cut separation from the first inhomogeneous magnon excitations 
(the absence of separation between the scaling laws would sign a quantum critical behaviour 



ii) the numbers of levels N$(N, S,p) in each 5* subset which for given N and S is determined 
by the number p of ferromagnetic sublattices of the underlying Neel order, 
in) the spatial symmetries of these N$(N, S,p) levels: The number and nature of the space 
group irreducible representations (S.G.I.R.) that appear in each S subspace are uniquely 
determined by the geometrical symmetries of the Neel order. 

The exact diagonalization of the Heisenberg hamiltonian on small samples enables us 
to examine these features and to determine the nature of the ordering. This approach has 
been used for the quantum Heisenberg antiferromagnet on the triangular lattice [f|,[7| and 
the J\ — J2 problem on the triangular lattice . 



III. THE LOW LYING LEVELS OF THE KAH 

The hamiltonian of the KAH, including further neighbor interactions, has the form: 

<i,j> «i,k» 

where the Sj are spin-1/2 quantum operators on the sites of the kagome lattice and < i, j > 
(resp. << i, k >>) denotes a sum over first (resp. second) neighbors. In this section, we shall 
only consider first neighbor exchange interactions (J2 = 0). At first sight, the spectrum of 
the KAH has one feature which is common to all the Heisenberg antiferromagnet spectra that 
we have been studying: the absolute ground-state has a total spin S = or 1/2 depending on 
the number of sites of the sample (N = 9, 12, 15, 18, 21, 24, 27) and the ground-state energy 
Eq(S) of the S sectors order with increasing S values. Lieb and Mattis [i0| have shown that 



this result is exact for bipartite lattices: our numerical results tend to indicate an extremely 
robust property (the theorem seems to be true for the Heisenberg antiferromagnet on the 
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triangular and kagome lattices and for the J\ — J2 problem: none of these situations can be 
reduced to a bipartite problem). Taken apart this feature the spectra of the KAH appear 
totally different from the spectra of the TAH. In Fig. 0, the low lying levels of the TAH and 
of the KAH are shown for the 27 sites samples: 

i) whereas the "Pisa Tower" is easily seen on the TAH spectrum, well separated from the 
magnons spectrum, there is absolutely no such scales in the KAH one, 

ii) the effective dynamics of the low lying levels of the KAH spectrum do not scale as 
S(S+1), 

Hi) the symmetries of the lowest lying levels of each S subspace do not allow the description 
of an ordered structure: for iV = 27, all the S.G.I.R., but one, appear in the low lying 
doublet states below the first S = 3/2 eigenstates, whereas their number N s and nature are 
strictly determined in the case of an ordered solution. 

One could argue that the proliferation of these low lying levels are the quantum counterpart 
of the infinite degeneracy of the classical ground state with respect to local spin rotations. 
The real question is: do the quantum fluctuations show any trend to select a specific Neel 
order? 

We have looked to this question for the so-called q = order, the \/3 x y/3 order and for 
any planar order (see section V). The q = order is studied in Ref. |4l| . We give here the 
details of the analysis concerning the \/3 x \/3 order, which is the favored solution found 
in the semi classical approaches |T8|-|22]|. The smallest lattices where periodic boundary 
conditions are compatible with this order are N = 9, 27 and 36 sites. In this section, we 
shall consider explicitly the N = 27 sample since the N = 9 sites is too small and the N = 36 
sample is too large to compute all the levels in each S sector. The QDJS associated with 
the a/3 x v3 state are homogeneous on each ferromagnetic sublattices (their wave-vectors 
are either k = or k = ±k : corners of the Brillouin zone). They do not break the C$ v 
symmetry of the lattice. The three irreducible representations characterizing the v3 x a/3 
order are the following: 
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T 1 :[k = 0, = n 27T/3 V = V, a x ^ = *] 

r 2 : [k = 0, K„V = ft 27r/3 # = = *] (3) 

^ r 3 : [k = ±ko, n 27T/3 v = #, o- a * = , 

where 7^ is a rotation of angle and er^ denotes an axial symmetry. The numbers 
Ns (N, S, p = 3) of levels in the "Pisa Tower" for each value of the total spin are given 
by the coupling of three N/6 spins: 

N/6 N/2 

piv/e ^ v n/q q piv/6 = (2S + 1) V s + (N/2-S + l)V s (4) 

S=0 S=N/6+l 

where T> s denotes the irreducible representation for a spin S. Therefore, one obtains the 
numbers N$(N, S,p = 3): 

N s (N,S,p = 3) = (2S + l)min(2S + l,N/2- S + 1). (5) 

We notice that in each S sector, an ordered solution countains a number of levels which is 
strictly related to the number of sublattices of the selected order: in the lower S subspace 
this number is independent of the sample size for p < 3 (it is the Hilbert space dimension of 
a rotator or a symmetric top for p = 2 (respectively p = 3)). In each S subspace, amongst 
the N S (N, S, p = 3) levels, the number of appearence (n^;) of the Tj IR can be computed 
following Refs. 0,^9]: 

4? = gZ)' & (^)Xi(*) Jv d(*) ( 6 ) 

k 

where the summation index k runs through the classes of the 5*3 group (isomorphic to C 3v ); 
Xi (k), N e \ (k) denotes respectively the characters of the IR and the number of elements 
in the class k (see Table |). The traces of the permutations of S3 in the S-subspace, denoted 
Tr (R k \ s ), are determined as: 

Tt(R k \s)=Tt[R k )-Tr(R k ), (7) 



s z =s 



S Z =S+1 



and 



Tr( (A,B) 
Tt((A,B,C) 



N/6 

Tr ( Li \ = ^ 5 t+v+x ,s z 

t,v,x=-N/6 
N/6 

t,v=-N/6 
N/6 

= ^ $3t,S z , 
t=-N/6 

where (A, B) (respectively (A,B,C)) stands for a two-body (respectively three-body) per- 

( s) 

mutation of S3. The final results of the computation of the rip are given for the N = 27 
sample in Table |I[ 

The lowest levels in the first 5* subspaces of the N = 27 sample spectrum are given in 
Table ffTTL The levels which have the good symmetries to describe a x \f?> antiferromagnet 
are displayed with an asterisk: most of these Ns(N = 27, S,p = 3) levels are rather far 
in the spectrum and many levels belonging to other S.G.I.R. proliferate between them. 
In fact, in the odd N samples, the number of low lying S = 1/2 states below the first 
S = 3/2 states increases very rapidly with the system size, the trend seems the same in 
the even samples (see Fig. |||). Altogether these numbers of low lying levels grow seemingly 
roughly as a N with a ~ 1.18 (resp. a ~ 1.14) in the odd (resp. even) samples. Note that 
these numbers lay between the ground-state degeneracy of the three states Potts model |42| 



(a ~ 1.134) and the degeneracy of the Dimer model f|3| (a = 2 1 / 3 ~ 1.26). This exponential 
proliferation of low lying levels with all spatial symmetries is certainly the deepest proof of 
the absence of long range antiferromagnetic order: it signs both the absence of a finite 
number of ordered sublattices (that is the absence of an antiferromagnetic order parameter) 
and the impossibility of a Neel symmetry breaking. 

The selection of order by quantum fluctuations previously observed in the J\ — J2 model 



on a triangular lattice is quite different |[29|| . In that last model the four-sublattice QDJS 
family is perfectly pure on the smallest non frustrating sample spectra: the two-sublattice 
QDJS family is a subset of the first family. Quantum fluctuations just stabilize this subset 
relatively to the entire four-sublattice family and thus build a simpler structure with an 



S 



order parameter of higher symmetry: They do not create the order parameter, but just 
renormalize it and increase its intrinsic symmetry. 

IV. ORDERING WITH A SECOND NEIGHBOR FERROMAGNETIC 

EXCHANGE COUPLING 

In order to ascertain our conclusion and reinforce the credibility of the method, we have 
studied with the same protocol the problem of the kagome lattice with a first neighbor an- 
tiferromagnetic interaction J\ and a second neighbor ferromagnetic interaction J 2 (see Eq. 
(§)) favouring the existence of a a/3 x v3 order. For large enough J 2 < 0, the spectrum 
has all the expected features of a "Pisa Tower" of QDJS (dynamics, number of states and 
symmetries) associated with this order (see Fig. ^ for the N = 9 sample). When \Jz\/Jx 
decreases the "Pisa Tower" disappears, indicating the existence of a quantum phase transi- 
tion. The estimate of the critical value of J 2 j J\ is a difficult task requiring a study of finite 
size effects. The smallest sizes compatible with the a/3 x v3 order are N = 9, 27, 36 sites: 
computing time and memory requirements for such an extensive study remain prohibitive. 
However, using appropriate twisted boundary conditions (±2-7r/3 around the z-axis, see Ref. 
and Appendix A), we can use the intermediate N = 12, 21 samples. The twisted bound- 
ary conditions break the rotational spin symmetry of the hamiltonian and fix the Neel plane 
perpendicular to the twist axis. The helicity of Neel order is thus fixed along the z-axis and 
the free dynamics of the system reduces to the precession of the total spin around z. The 
effective hamiltonian reads: 

H eS = (9) 

where 1% is the inertia of the top along the z-axis. The degeneracy of each Si subspace is 
2 (±S Z ). The IR characterizing the a/3 x a/3 order is Ti : [k = 0, 7^/3* = ^ = *] . 
Fig. |5| shows the tower of states associated with the a/3 x a/3 order for the iV = 21 sample 
with | J 2 / J\\ = 1. 
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The order parameter or the spin stiffness of these levels scale as iV -1 / 2 , whereas the slope 
of the "Pisa Tower" scales as N~ 2 . At first sight, it seems thus more efficient to use this 
information to look at the transition between order and disorder as \ Ji\/Ji is decreased. 
This can be done in a rather naive way by looking at the spin gap between the ground-state 
of ^min and S Z)Tain + 1: in the ordered phase this quantity should scale as iV _1 (respectively 
N~ l l 2 and N° in the critical and disordered phase): this study is done in Fig. [|for the case 
I Jij J\ I = (pure K AH) . It shows thanks to the trends in the larger sample sizes that the 
KAH is certainly not ordered, probably not critical but truly disordered. 

This naive use of the "Pisa Tower" does not account for the whole qualitative information 
contained in the QDJS family which is described both by the effective hamiltonian and by 
the space group symmetry of the levels. This information is incorporated in the index R 
measuring the "degree of order" present in the low lying levels, and defined as follows: 

i) First, we consider the ordered phase (with \J^/ J\\ > 1) and concentrate on all the low 
lying levels of the "Pisa Tower" that are lower in energy than the softest magnons: for each 
sample size this determinates the number of S z sectors which give a consistent picture of the 
"quasi-classical ground-state". Precisely we search for the lower level of the spectrum with 
a non zero wave vector in the magnetic Brillouin zone (£' m j n (k)) and determine S z ^ max as 
the largest value of S z such that -^(^max + 1) < E min (k). By definition E (S z>max + 1) is 
thus smaller than the energy of the softest magnon of the sample and in the thermodynamic 
limit S^max grows roughly as iV 1 / 4 . 

ii) In each S z subspace considering all the levels in the energy range [E (S Z ),E (S Z + 1)] 
we define r$ z as the ratio of the number of levels compatible with y/3 x y/3 order to the 
total number of levels of this range (r$ z is a number between and 1). 

in) We then compute the index R measuring the degree of order of the sample as the 
average of the ratios r Sz for S z running from (resp. 1/2) up to S^max- 

This index is equal to one if the lowest levels form a true "Pisa Tower" and decreases 
when other IRs appear in the low lying levels of the spectrum when \J2/Ji\ is decreased. R 
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is thus a measure of the breakdown of order which includes qualitative information on the 
low lying levels. It stands on the quantities that scale the more rapidly with the system size. 

The variations of this index as a function of \J2\fJ1 and of the sample size are given in 
Fig. [7]. When \Jz/ J\\ = this index goes very rapidly to zero with N: this conforts the 
idea that the KAH is indeed disordered and that quantum fluctuations show no tendency 
to select "order from disorder" in the disordered phase. More unexpected, the finite size 
scaling on R indicates that a small ferromagnetic exchange is not sufficient to establish long 
range order in the system and that the value of the critical ratio (| Ji\/ J\) c is probably larger 
than 0.5. 

V. INCOMMENSURATE MAGNETIC ORDER 

The previous study discards the hypothesis of a \f?> x y/3 order in the pure KAH (we did 
the same check with the same conclusion for the q = order in Ref. pifl). Using twisted 
boundary conditions S ri+Tct=1 2 = 7^ z ($ Q , = i i2 )S ri across the sample defined by the vectors 
T a =i.2 (see the appendix for more details), we have extended our search to any planar 
antiferromagnetically ordered configurations, either commensurate or incommensurate. The 
existence of a planar order, if any, would be signed by a minimum of the ground-state energy 
for a given pair of twist angles ($1, $2) and the appearence of a "Pisa Tower" for this couple 
of parameters [|7| . A typical result of a set of diagonalizations is shown in Fig. [8] for the 
N = 21 sample. Sweeping the Brillouin zone for ($1, $2), we have studied in this way the 
spectra of the N = 9, 12, 15, 18, 21, 24, 27 samples. We observe the following properties: 

i) the influence of the twisted boundary conditions is very small, much smaller than for the 
TAH: on the N = 21 sample the effect of the twist on the ground-state of the KAH is only 
8% of the same effect on the TAH. This is coherent with the picture of a disordered, liquid 
system, 

ii) we do not find a signature of any planar antiferromagnetic order for any size. Very 
shallow minima appear in the spectra, but they are never associated with a tower of QDJS 
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and the position of these minima changes from place to place with the sample size, 
in) the "spin-gap" between the ground-state energy of the S z = (or 1/2) subspace and 
the ground-state energy of the S z + 1 subspace (AE S = E (S Z + 1) — E (S Z )) has only small 
variations with the twists (Fig. |9[). These variations appear systematic, and show a different 
trend in the even and odd samples: this could be related to the fact that the odd samples 
only accomodate spin- 1/2 excitations of the thermodynamic absolute ground-state which is 
a true singlet. (This hypothesis is examined in a companion paper p4|). 



VI. CONCLUSION 

Using the analysis of the low lying levels of the Heisenberg antiferromagnet on a kagome 
lattice we have shown new evidences that the system has no planar antiferromagnetic long 
range order at T = 0. Introducing a small second neighbor ferromagnetic exchange coupling 
does not seem to be sufficient to establish long range order: from the experimental point of 
view this is good news as it could perhaps enlarge the number of candidates for a spin liquid 
behaviour. 

The theoretical study of the tower of low lying levels ( "Pisa Tower" ) that we have devel- 
opped in this paper seems potentially useful to give an approximate location of the transition 
from order to disorder even on small samples: its advantage on other approaches stands on 
the finite size scaling of the parameter we are looking at and on the inclusion in this param- 
eter of qualitative information on the macroscopic ground-state. 

According to our present results, this spin-1/2 model exhibits a quantum critical point 
at a non zero Ity/Ji- It would be interesting to investigate the universality class of this 
quantum critical point and see how it may compare to the theoretical predictions of the non 



linear sigma model for canted antiferromagnets j45 . 
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APPENDIX A: TWISTED BOUNDARY CONDITIONS 



For a sample defined by the two vectors: 

Ti = 2(1 + m)ui + 2mu 2 

(Al) 

T 2 = 2/ux + 2(1 + m)u 2 

where ui and u 2 are two unit vectors of the kagome lattice, and m and n are two integers 
related to the number of sites of the sample by: N = 3(l 2 +m 2 +lm), the boundary conditions 
are defined through: 

Sp i+ T a=1 , 2 = ft,($a=i l2 )S Pi (A2) 

where denotes the sites of the kagome lattice. In order to recover the translation invariance 
that seems to be broken by these boundary conditions, the spin frame at the point + ui 
(resp. Yi + u 2 ) is rotated with respect to the spin frame at point by an angle Q\ (resp. 
6* 2 ). The boundary angles $ a =i, 2 are related to 6 a =\,2 by the relations: 

$! = 2(l + m)e 1 + 2m6 2 

(A3) 

$ 2 = 2W 1 + 2(1 + m)6 2 . 
The hamiltonian in the new frame reads: 

H = 2.h J2 S ri .^(^)S ri+UM , (A4) 

i=l, N 

/i=l,3 

Q\ and 9 2 are changed step by step so that $i and $ 2 sweep the appropriate fraction of the 
"Brillouin zone" of this problem . 
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TABLES 
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TABLE I. Character table of the permutation group S3. First line indicates classes of permu- 
tations. The number of elements in each class is iV e j. 
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TABLE II. Number of occurrences n v . of each irreducible representation Lj (i = 1,2,3) with 
respect to the total spin S. 
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TABLE III. Lowest eigenenergies (with degeneracy and quantum numbers) of the N = 27 
sample spectrum. Components of k are in units 6ir/N. In the 3 last columns, 1 stands for 
invariant under the symmetry, for no symmetry and -1 for a non trivial phase factor under the 
above mentionned symmetry (i.e. e* 27r / 3 for the rotation of 27r/3 and -1 for the two others). Stars 
indicate states possessing the symmetries associated with the \/3x\/3 state. The double horizontal 
bar indicates the omission of 127 S = 1/2 states before the first levels in the S = 3/2 subspace. 
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FIGURES 

FIG. 1. Two classical planar states of the KAH: the \/3 x a/3 and q = states. 

FIG. 2. The low lying energy levels of the TAH and KAH spectrum of the N = 27 sample. 
The levels which possess the symmetry expected for an ordered solution are denoted by a star. 
The "Pisa Tower" of the TAH is easily seen, well distinct from the first magnon excitations. In 
the KAH on the contrary the levels candidate for the building of a tower of states are mixed with 
other representations in a continuum of excitations. 

FIG. 3. Logarithm of the number An of S = 1/2 levels below the first S = 3/2 level as a 
function of the sample size iV (black triangles). This number which does not take into account the 
two-fold magnetic degeneracy is compared to the same quantity (square symbols) for the even N 
samples (i.e. number of S = levels below the first S = 1 level ). 

FIG. 4. Low lying spectrum of the N = 9 sample with \ Ji\l J\ = 1 and periodic boundary 
conditions, as a function of S(S + 1). Note the "Pisa Tower" associated with the \/3 x \/3 state. 

FIG. 5. Low lying spectrum of the iV = 21 sample with | Ji\l J\ = 1 as a function of S 2 . A twist 
of 27r/3 in the boundary conditions is applied to accomodate the y/3 x \/3 state with the sample 
size. Due to the boundary conditions, the "Pisa Tower" reduces to one IR for each S z value. 

FIG. 6. Finite size study of the "spin gap" AE S = Eo(S Z)m i n + l) — Eo(S Zym in) of the pure KAH 
(J2 = 0) as a function of the sample size. This energy is 0(N a ) with a = —1, —1/2 or whether 
the system is ordered, critical or disordered. Continuous lines (resp. broken lines) are guides for 
the eye through the even N (resp. odd N) results. The comparison of the three quantities AE S , 
N 1 / 2 x AE S and N x AE S versus iV 1 favors the hypothesis of spin disorder. 

FIG. 7. Behavior of the "index of order" (R) as a function of | J2I/J1 and of the system size. 
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FIG. 8. Variation of the energy per link (2Sj.Sj) of the low lying levels of the N = 21 sample 
versus twisted boundary conditions <E>i ($2=0). Oj are the points where 3>i= (mod 2tt). O-j is 
the first point where the twist per link is (mod 2n). The points (3>i=7r (mod 2tt)) are in the 
middle of Oi and Oj+i- Because the figure is symmetric with respect to ^3, the part A3 — O7 is not 
represented here. Note that most of the levels do not come back to their original assignation after 
a 2tt twist of the boundary conditions. In fact on this small size sample, due to an extra-symmetry 
when <I>i=0, the uniform k=0 ground-state is degenerate with the first star of k / eigenstates. 
Only the k=0 states and their continuation have been shown in OA\. The 6 other k / in OA 
are found by folding the 0{A{ onto OA. Full lines stand for states going continuously to k=0 
chiral states (complex IRs of C3), dashed lines stand for states going continuously to the first star 
of k / eigenstates the k=0 non chiral states (trivial IR of C3); bold line: the first S z =3/2 level. 

FIG. 9. Variations of the "spin gap" with boundary conditions: the small horizontal tick gives 
the value of the gap for periodic boundary conditions. There is a systematic effect: the AS" = 1 
excitation energy decreases with the twists of the boundary conditions in the even N samples and 
increases in the odd N ones. 
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